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ORIGINAL ARTICLE 

Application of a Physiologically Based Pharmacokinetic 
Model to Predict 0ATP1B1 -Related Variability in 
Pharmacodynamics of Rosuvastatin 

RH Rose\ S Neuhoff\ K AbduljaliP, M Chetty\ A Rostami-Hodjegan^'^ and M Jamei^ 

Typically, pharmacokinetic-pharmacodynamic (PK/PD) models use plasma concentration as the input that drives the PD model. 
However, interindividual variability in uptake transporter activity can lead to variable drug concentrations in plasma without 
discernible impact on the effect site organ concentration. A physiologically based PK/PD model for rosuvastatin was developed 
that linked the predicted liver concentration to the PD response model. The model was then applied to predict the effect of 
genotype-dependent uptake by the organic anion-transporting polypeptide 1 B1 (0ATP1 B1) transporter on the pharmacological 
response. The area under the plasma concentration-time curve (AUCg J was increased by 63 and 111% for the C.521TC and 
C.521CC genotypes vs. the C.521TT genotype, while the PD response remained relatively unchanged (3.1 and 5.8% reduction). 
Using local concentration at the effect site to drive the PD response enabled us to explain the observed disconnect between the 
effect of the 0ATP1B1 c521T>C polymorphism on rosuvastatin plasma concentration and the cholesterol synthesis response. 
CPT Pharmacometrics Syst. Pharmacol. (2014) 3, e124; doi:10.1038/psp.2014.24; published online 9 July 2014 



Physiologically based pharmacokinetic (PBPK) models are 
increasingly being used to predict the impact of physiological 
and pathophysiological patient factors and concomitant med- 
ication on drug exposure to support drug development and 
regulatory submissions.^'^ Typically, the end point of a PBPK 
model is a prediction of the pharmacokinetics of a drug, while 
the ultimate success of a drug is dependent on demonstra- 
tion of efficacy without toxicity. Because PBPK models can 
predict drug concentration in tissues and in plasma, a natural 
progression is to link them to pharmacodynamic (PD) mod- 
els via the concentration at the site of action. ^'^ Compared 
with the traditional approach of pharmacokinetic/pharmaco- 
dynamic (PK/PD) modeling that uses plasma concentration 
to drive the response, this may allow a better understand- 
ing of true PD variability vs. variability that results from drug 
disposition to the site of action.'' This is particularly pertinent 
where transporters are involved in drug disposition to its 
effect site. In this case, interindividual variability in transporter 
activity can result in a lack of correlation between plasma 
concentration and concentration at the site of action between 
individuals.^ 

Rosuvastatin is a 3-hydroxy-3-methylglutaryl coenzyme 
A reductase inhibitor, thus reduces the conversion of acetyl- 
CoA to mevalonic acid (MVA), which is the rate-limiting 
step in hepatic cholesterol biosynthesis.*^ Rosuvastatin has 
low passive permeability across biological membranes that 
limits distribution to tissues and oral absorption.^ However, 
rosuvastatin is extensively distributed into the liver, its major 
site of action, through the action of specific uptake transport- 
ers, including the organic anion-transporting polypeptides 
0ATP1 B1 , 0ATP1 B3, and 0ATP2B1 and the sodium-depen- 
dent taurocholate cotransporting polypeptide. ''■^ Both liver 
canalicular and intestinal efflux of rosuvastatin are medi- 
ated by the breast cancer resistance protein.'" Multidrug 



resistance-associated protein-2 also contributes to the liver 
canalicular efflux of rosuvastatin but plays a more significant 
role in rats than in humans.'' '^ The liver is a site of elimina- 
tion of rosuvastatin, predominantly through biliary elimination 
and to a lesser extent metabolic elimination (-10%).'^ 

Genetic variants of 0ATP1B1 and breast cancer resis- 
tance protein have been identified that contribute to interin- 
dividual variability in rosuvastatin disposition, exposure, and 
therapeutic or side effects."'^ In this study, we focus on the 
0ATP1B1 c.521T>C single-nucleotide polymorphism (SNP). 
This SNP has been associated with increased exposure to 
rosuvastatin because of reduced clearance'"'^ and a sig- 
nificant increase in the risk of myopathy for the statins sim- 
vastatin and atorvastatin, although a statistically significant 
increase in risk has not been found for rosuvastatin.'**"^" 

A PK/PD model describing the effect of rosuvastatin on 
plasma MVA concentration has been published that uses an 
indirect response model with a circadian rhythm on the input 
rate.^' This model is typical of PK/PD models in that it uses 
the total plasma concentration to drive the PD model. How- 
ever, the concentration of rosuvastatin at the site of action, 
i.e., the hepatic unbound intracellular water concentration 
(CU|^,,), is a more relevant driving concentration for the PD 
model. This is supported by a recent publication that showed 
an improved correlation in the cholesterol-lowering effect 
between humans and a mouse model when hepatic extrac- 
tion was accounted for.^^ PBPK models have previously been 
described for rosuvastatin that account for transporter-medi- 
ated disposition and allow prediction of hepatic Cu,^."'^^'^" 

The aim of this study was to demonstrate the added util- 
ity of linking the PBPK model-predicted concentration at the 
site of action to the PD response compared with the plasma 
concentration, using rosuvastatin as an example. PBPK and 
PD models were integrated, creating a PBPK/PD model that 
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Figure 1 Schematic of the physiologically based pharmacokinetic-pharmacodynamic (PBPK/PD) model used in this study. The unbound 
concentration of rosuvastatin in the intracellular water of hepatocytes (CU|„) and the impact of liver uptake transporter activity are predicted 
by a permeability-limited liver model within a full PBPK model. This accounts for the distribution of unbound drug from the blood into the 
extracellular water and the passive diffusion and active transport {K,|;j^ ,^ and K,|„ J across the sinusoidal membrane between the intracellular 
and extracellular space. The impact of ionization is considered (red circles represent ionized drug and blue circles unionized drug) as well 
as the association of drug with extracellular proteins and intracellular acidic phospholipids and partitioning into neutral lipids and neutral 
phospholipids. Elimination of intracellular drug via biliary clearance and metabolism are considered. The PD response model uses a modified 
indirect response model input, with the drug effect resulting from inhibition of the input rate (KJ and driven by the predicted liver CU|^^. 



links the unbound hepatocellular concentration, predicted by a 
permeability-limited liver model within a full PBPK model, to the 
rate of cholesterol synthesis over time (Figure 1). The devel- 
oped model was then used to predict the impact of 0ATP1B1 
C.521TT, TC, and CC genotypes on the PK and PD of rosuvas- 
tatin, and to compare the predictions with those when plasma 
concentration was used to drive the PD response and with clini- 
cal data. Simulations also investigated the impact of reduced 
hepatic uptake transporter function on plasma, liver, and mus- 
cle concentration, which relates to the myopathy side effect of 
statins, in addition to the rate of cholesterol synthesis. 

RESULTS 

The 0ATP1B1 hepatic uptake clearance (CL^,.^) for the 
0ATP1B1 C.521TT, TC, and CC genotypes was fit using 
published plasma concentration-time data for subjects 
stratified by this genotype to obtain values of 126, 30, and 
0 ^il/min/10'' cells, respectively. Parameter estimates were 
able to satisfactorily recover the clinical plasma concentra- 
tion profiles for subjects grouped by genotype (Figure 2). 



For all three genotype groups, the mean observed maxi- 
mum concentration (C^^J and area under the plasma 
concentration-time curve (AUC(,_^g,,) were within the range 
of those for the 10 simulated profiles matched in terms 
of subject age, proportion of females, and study size for 
each genotype (Supplementary Table SI). However, in all 
cases, the simulated time after administration to maximum 
plasma concentration (7"^^,^ underestimated the observed 
data (Supplementary Table S1). The simulations predict 
that the mean AUC„_„ and C^^^ were increased by 86 and 
90% for the C.521CC genotype and 62 and 60% for the 
C.521TC genotype, relative to the C.521TT genotype. The 
clinically observed increases in AUC„ and C were 

^ 0-«' max 

reported as 62 and 79%, respectively, for the c.521 CC gen- 
otype group and 57 and 52%, respectively, for the c.521 TC 
genotype group. 

The estimated ICj.^ and Hill coefficient for the drug effect 
based on rosuvastatin liver Cu^^ were 0.13 jjmol/l and 1.3, 
respectively. The final model incorporating fitted the param- 
eters allowed adequate recovery of PK and PD profiles 
(Figure 3). The simulated reduction in MVA AUC relative to 



CRT: Pharmacometrics & Systems Pharmacology 



PBPK Model of OATPIBI-Related PD Variability 

Rose efal. 




Time (h) Time (h) Time (h) 



Figure 2 Simulated and observed plasma rosuvastatin concentration profiles for the (a) wild-type (c.521TT), (b) heterozygous, and (c) 
homozygous deficient (c.521 CC) organic anion-transporting polypeptide 1 B1 (0ATP1 B1 ) genotypes. Gray lines represent simulated individual 
trials, and the solid black line represents the simulated mean of 10 trials. Circles represent data extracted from Pasanen ef a/..^^The simulation 
study design was matched to that of the clinical study. 




Figure 3 (a) Simulated and observed plasma rosuvastatin concentration profile after multiple daily dosing (10mg) for 14 days. Gray lines 
represent simulated individual trials, the solid black line represents the simulated mean of 10 trials, and circles are mean observed data, (b) 
Mean simulated (lines) and observed (markers) plasma MVA profile before rosuvastatin administration (baseline) and after morning (07:00) 
or evening (18:00) rosuvastatin administration. The time is 06:00 at 31 2 h. Observed data are from IVIartin ef aA.^The simulation study design 
was matched to that of the clinical study. MVA, mevalonic acid. 



baseline was 34 and 27%, respectively, for the evening and 
morning dose of rosuvastatin. 

Simulations were performed to predict the effect of the 
0ATP1B1 c.521T>C sequence variation on the PD response 
using the established PBPK/PD model (Figure 4). Analysis 
of the simulated data showed that the mean rosuvastatin 
plasma AUCq__ was increased by 63 and 111% for the het- 
erozygote and CC-homozygote genotypes relative to the wild 
type (Table 1).The liver Cu^^ AUC„_„ was reduced by 5.7 and 
9.6%, respectively (Table 1).The average MVA AUG relative 
to baseline corresponding to these sequence variations was 
increased by 30 and 35% when plasma concentration was 
used as the input to the PD response model, as in the original 
modeP^ (Table 1). However, when liver Cu^^ was used as the 
input to the PD response in our modified model, the simu- 
lated average MVA AUG relative to baseline was reduced 
by 3.1 and 5.8%, respectively (Table 1). The latter model 
is more consistent with clinical data and with mechanistic 
understanding of statin action (see Discussion), so this was 
used in further analysis. 

There was a large interindividual variability in both the 
plasma AUGq^ of rosuvastatin and the reduction in MVA AUG 



from baseline (Figure 5), resulting in overlap between the 
interquartile ranges for the three 0ATP1 B1 genotype groups. 
This overlap is most pronounced for the PD response, in 
which median values for each group are within the interquar- 
tile range for all groups. 

Further analysis of the impact of the total hepatic uptake 
transporter GL^^.^ on the plasma, muscle, and liver Gu^^ rosu- 
vastatin AUGq_24^ and MVA AUG relative to baseline predicted 
by the integrated PBPK/PD model was performed using sen- 
sitivity analysis for a range between 0 and 250 ^jl/min/1 0** cells. 
Both plasma and muscle exposure to rosuvastatin decreased 
as the overall hepatic uptake GL^^^ increased (Figure 6a, c). 
In contrast, both liver Gu,^ and PD response increased as the 
overall hepatic uptake CL^,.^ increased (Figure 6b, d). In all 
cases, sensitivity to the value of CL^^J is greatest when the 
value of CL^^J is low. The elasticity index (El) (Figure 6e) is 
a measure of the relative change in selected variable to the 
relative change in the GL^^.^, normalizing the scale for com- 
parison of the sensitivity of multiple model output parameters. 
The El of the muscle rosuvastatin AUG to uptake transporter 
GL.^,.^ is identical to the plasma AUG and has its greatest 
absolute value, indicating the greatest proportional change 
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Figure 4 The effect of the organic anion-transporting polypeptide 1 B1 (OATP1 B1 ) sequence variation on the simulated mean (a) plasma 
concentration and (b) liver unbound intracellular water concentration (CUn^) of rosuvastatin and (c) plasma mevalonic acid (MVA) concentration 
using plasma concentration or (d) liver Cu^^ as the driving concentration for the response. Data are the mean of 100 simulated individuals 
based on a population that was 50% female with an age range 20-50 years. Individuals were dosed with 10mg oral rosuvastatin at 18:00 
daily for 5 days. The time is 18:00 at 108h. 



Table 1 Mean simulated plasma and liver exposure and pharmacodynamic response to rosuvastatin for the three OATP1 B1 genotypes investigated 



0ATP1B1 genotype 


Plasma AUC„__ 
(ng/ml-h) 


Liver Cu,„AUC^ 
(ng/ml-h) 


MVA 24 h AUC relative to baseline 
(%): plasma concentration input 


MVA 24 h AUC relative to 
baseline (%): liver Cu,^ input 


C.521TT 


35 


120 


35 


36 


C.521TC 


57 


114 


46 


35 


C.521CC 


74 


109 


51 


34 



Mean simulated plasma area under the curve (AUC„ J, liver unbound concentration in intracellular water (Cu^^) AUCj, and mevalonic acid (MVA) AUC relative 
to baseline of rosuvastatin for the three organic anion transporting polypeptide 1 81 (OATP1 B1 ) genotypes investigated. MVA AUC was simulated by using either 
plasma concentration or liver Cu,^ as the input to the PD model. Data are the mean of 100 simulated individuals based on a population that was 50% female with 
an age range 20-50 years. Individuals were dosed with lOmg oral rosuvastatin, either at single dose {plasma and liver CU|„ AUC,^) or at 18:00 daily for 5 days 
(MVA AUC, results for final dose). 



in AUC, when hepatic uptake transporter CL^,.^ is between 
175 and 225 jjl/min/10^ cells. Changes are in the opposite 
direction for liver CU|^,, and PD response (positive rather than 
negative), with greatest elasticity at overall hepatic CLi^,.^ of 
25 ^il/min/10^ cells. The PD response shows a lower elasticity 
to changes in CL^,.^ than the liver Cu,^. 

DISCUSSION 

This study aimed to demonstrate the added utility of link- 
ing the PBPK model-predicted concentration at the site of 



action to the PD response compared with that of the use 
of plasma concentration. Rosuvastatin was selected as an 
example drug because there are clinical data demonstrating 
the contrasting effect of hepatic transporter activity, specifi- 
cally for the 0ATP1B1 c.521T>C genotype, on the plasma 
concentration and PD response. In addition, both a PBPK 
model that includes hepatic transporter activity and a PD 
model for rosuvastatin have been published. Thus, using 
these models, a combined PBPK/PD model could be gener- 
ated with changes to only three parameter values to account 
for genotype-specific 0ATP1B1 uptake activity and altered 
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Figure 5 Comparison of the interindividual variability of simulated 
(a) plasma rosuvastatin AUCg _ and (b) reduction in MVA AUC 
from baseline for 0ATP1 B1 C.521TT, TC, and CC genotype groups. 
Box and whisker plots: horizontal lines from bottom to top are the 
minimum, 25th percentile, median, 75th percentile, and maximum 
values. The diamond symbol is the mean. Distributions obtained 
for 100 simulated healthy volunteer individuals for each genotype. 
Variability was included in PBPK model parameters but not the PD 
model parameters. 

sensitivity of the PD response to the different input concen- 
tration used. Compared with using plasma concentration as 
the driving concentration for the PD response, the model 
is better able to capture the clinical effect of the 0ATP1B1 
c.521T>C SNP on the therapeutic effect of rosuvastatin. 

The genotype-specific 0ATP1B1 CL^,.^ for rosuvastatin, 
considering only the c.521T>C SNP, was estimated using 
clinical data from Pasanen et al..^'' The fitted values of 
126, 30, and 0 Ml/min/10« cells for the C.521TT C.521TC, 
and C.521CC genotypes, respectively, are consistent with 
an average 0ATP1B1 CL^,.^ of 109 tJl/min/10'* cells deter- 
mined for rosuvastatin for the north European Caucasian 
population, in which the C.521T allele predominates. 
Simulations using the parameter estimates for the differ- 
ent 0ATP1B1 genotype groups and a simulated study 
design matched to that reported by Pasanen et al.''^ were 
able to recover the clinical data well (Figure 2). The simu- 
lated average increase in plasma AUC for the 0ATP1B1 
C.521CC genotype relative to the TT genotype (111%) is 
also in close agreement with that reported in a separate 
clinical study for white subjects, in which a 117% increase 
in activity was observed." However, estimation of complete 
loss of transport activity for the 0ATP1B1 C.521CC geno- 
type conflicts with in vitro studies that have shown reduced, 
but not complete loss of, rosuvastatin transport activity for 
expressed 0ATP1B1 with the C.521C SNP''=^'' Accounting 
for residual 0ATP1B1 activity for the C.521CC genotype 
would reduce the simulated plasma rosuvastatin concen- 
tration, which tends to overestimate mean observed data 
with current settings (Figure 2). 

Although the estimated CL„,^ for the 0ATP1B1 C.521TC 
genotype group was able to recover clinical data from Pas- 
anen et al." a much smaller increase in plasma AUC of only 
6.3% for the C.521TC relative to the C.521TT genotype was 
reported in another clinical study of white subjects." Thus, 
this study was poorly predicted by our model and the impact 



of the C.521TC genotype on plasma concentration is con- 
troversial. A large interindividual variability in rosuvastatin 
exposure is predicted for the different genotype groups 
(Figure 5), indicating the influence of many covariates, most 
of which remain unknown in clinical studies. Discrepancies 
between the impact of 0ATP1B1 c.521T>C sequence varia- 
tion between clinical studies and between clinical and in vitro 
data may reflect limitations of the fitting approach, which 
fixed parameter values for all but the fitted parameter and 
used average clinical data from a small study. The availability 
of more clinical plasma concentration data stratified by geno- 
type or data that are required for in vitro-in vivo extrapolation 
(such as absolute transporter abundance and in vitro activity 
data with extrapolation factors) would help generate CL^,.^ 
values with increased confidence. Recent reports on mea- 
surement of transporter protein abundance data in human 
hepatocytes can improve the in vitro-in vivo extrapolation of 
transporters.^^ 

The PD model was an adaptation of that reported by 
Aoyama et al." which describes the change in plasma MVA 
concentration in healthy subjects receiving lOmg oral rosu- 
vastatin daily. The model was modified by refitting the drug 
effect parameters (the IC^^ and Hill coefficient) when the 
hepatic Cu^^ predicted from the PBPK model was used as 
the input to the PD model. The estimated ICj.„ for inhibition 
of 3-hydroxy-3-methylglutaryl coenzyme A reductase activity 
by rosuvastatin (0.13 ^jmol/l) is considerably higher than the 
reported IC^^ for the inhibition of 3-hydroxy-3-methylglutaryl 
coenzyme A reductase activity of a purified catalytic frag- 
ment (3.5-5.4x10"^ |jmol/l) and in isolated rat hepatocytes 
(1.3x10"^ ^imol/l)."'" This suggests that the PD model may 
lack sufficient mechanistic detail to make in vitro-in vivo 
extrapolation of the IC^-g appropriate or that the in vitro meth- 
odology fails to reflect the in vivo situation. 

The final PBPK/PD model maintained the ability to describe 
the plasma MVA profile adequately (Figure 3). The simulated 
reduction in MVA AUC relative to baseline of 34 and 27%, 
respectively, for the evening and morning dose of rosuvas- 
tatin is in reasonable agreement with the reduction of 40 and 
32% simulated by Aoyama et al}^ using their PK/PD model 
with plasma concentration as the PD model input. The study 
used in fitting of the moder is the only published clinical study 
that has reported the change in plasma MVA concentration 
after rosuvastatin administration. Consequently, the PD 
model has not been tested against other datasets, including 
those using a dose other than lOmg rosuvastatin or in differ- 
ent ethnic or disease populations. Therefore, caution should 
be taken in generalizing the model to other populations or 
dosing regimens, particularly in terms of the precise quantita- 
tive description of the MVA profile. The model also assumes 
that there is no influence of 0ATP1 B1 activity on the concen- 
tration of MVA at baseline. An increased baseline cholesterol 
synthesis rate has been reported for the 0ATP1 B1 C.521CC 
genotype, although there was no effect on the total plasma 
cholesterol.^" Mechanistically, the increased cholesterol 
synthesis rate may be related to reduced OATP-mediated 
hepatic bile acid uptake, leading to reduced hepatic bile acid 
concentration and removal of the inhibitory effect of hepatic 
bile acid on cholesterol catabolism.^" However, data are for 
a small study size (32 subjects, only 4 with the C.521CC 
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Figure 6 Sensitivity analysis of the influence of total uptake transporter intrinsic clearance (CL^,^) on (a) plasma area under the curve 
(AUC„ 2^,^), (b) liver unbound concentration in intracellular water (Cu^^) A\JC^^^^, (c) muscle AUCj, of rosuvastatin, and (d) the reduction 
in plasma mevalonic acid (IVIVA) AUG relative to baseline, (e) The elasticity index allows direct comparison of the relative change of plasma, 
liver CU|^, and muscle AUG of rosuvastatin and reduction in MVA AUG ratio to the relative change in total uptake transporter GL^,^ MVA, 
mevalonic acid. 



genotype), and the biomarker for cholesterol synthesis was 
not the plasma MVA concentration, so this was not included 
in the present model. 

When rosuvastatin liver CU|^,, was linked to the PD response, 
the PBPK/PD model predicted a small reduction in the MVA 
AUG relative to baseline for individuals with the 0ATP1B1 
c.521T>C SNP (3.1% for C.521TC and 5.8% for C.521CC). In 
contrast, a larger increase in the MVA AUG relative to base- 
line (30% for C.521TC and 35% for C.521CC) was predicted 
if the PD response was linked to the plasma concentration. 
The former is in agreement with a number of studies that 
have shown either no effect or a slightly reduced therapeutic 
response associated with the 0ATP1B1 c.521T>C SNP for 
the cholesterol-lowering response to statins. ^''^"^^ For exam- 
ple, in a genome-wide study of over 3,000 patients allocated 
to receive rosuvastatin 20 mg daily for a year, the 0ATP1B1 
c.521T>C SNP was associated with a 2.6% reduction in 
fractional low-density lipoprotein cholesterol reduction per 
allele. Large interindividual variability in the PD response is 
predicted by the PBPK/PD model with considerable overlap 
in the distributions (Figure 5), suggesting that a large study 
size is required to detect the small effect of the 0ATP1B1 
genotype on the PD response to rosuvastatin. A limitation of 
this study is that interindividual variability was not included for 
the parameters describing the PD response because the pub- 
lished PD model adapted for this study was estimated from 
average profiles. ^'^^ As a result, variability is only introduced 
by the PBPK model parameters, and PD variability is likely 
to be underestimated, possibly to a large extent because PD 
variability can be much greater than PK variability. 



Taken together, the results indicate that reduced 0ATP1 B1 
activity results in a relatively large increase in plasma rosuv- 
astatin concentration and a small decrease in PD response. 
These results are in agreement with the predictions of Wata- 
nabe et al.^ for the effect of uptake transporter activity on 
plasma and liver concentration of pravastatin. Reduced 
uptake activity may be expected to decrease liver Cu^^, but 
this effect is countered by increased plasma concentration 
of rosuvastatin because reduced liver exposure leads to 
reduced hepatic drug elimination. The higher plasma con- 
centration results in increased liver unbound concentra- 
tion in extracellular water (Cu^^^), which drives both passive 
and active uptake into the liver. Because the PD response 
is driven by the liver Cu^^, both show a similar sensitivity to 
uptake transporter activity, but the relative sensitivity of the 
PD response is slightly lower because of the nonlinearity in 
the PD model (Figure 6e). A limitation of the study by Wata- 
nabe et al.^ is that predictions were based on sensitivity 
analysis to probe model behavior that was not confirmed by 
clinical or preclinical data. Our study provides verification that 
the modeling approach is able to recover clinical data for the 
impact of hepatic uptake transporter activity on the pharma- 
cokinetics (plasma concentration) and pharmacodynamics of 
rosuvastatin for specific 0ATP1 B1 genotypes. 

The El, a normalized measure of sensitivity, for the effect of 
uptake transporter activity on muscle AUG exactly matched 
with that of the plasma AUG (Figure 6e). This is expected 
because prediction of muscle concentration was based on 
the perfusion-limited model, thus assuming uptake by rapid, 
passive diffusion into the tissue. The results are in agreement 



CRT: Pharmacometrics & Systems Pharmacology 



PBPK Model of OATPIBI-Related PD Variability 

Rose efal. 



with the association that has been observed between the 
plasma concentration of statins and the risl< of muscle-related 
side effects, such as myopathy and rhabdomyolysis.^" Thus, 
in contrast to the results for the cholesterol-lowering effect 
of rosuvastatin, this supports the use of plasma concentra- 
tion as a surrogate for the concentration at the site of action 
(muscle) in assessing risk of statin-induced muscle toxicity 
in individuals with different hepatic uptake transporter activ- 
ity. One study has suggested a role for transporter-mediated 
uptake, by 0ATP2B1, and efflux in skeletal muscle expo- 
sure to rosuvastatin.^^ However, at present, the importance 
of transporter-mediated uptake of rosuvastatin into muscle 
remains unclear. Expression of mRNA for 0ATP2B1 was 
considerably lower in skeletal muscle than in the liven^*^ and 
insufficient data are available to model the impact of specific 
transporters on rosuvastatin uptake into muscle. 

The potential power of linking PBPK to PD models has pre- 
viously been recognized;^ however, there are few published 
examples that demonstrate successful application of this 
approach and the added value it can offer. To our knowledge, 
this is the first published study to use the liver concentra- 
tion from a full PBPK model to drive the PD response and 
demonstrate improved ability to assess the impact of trans- 
porters involved in the uptake to the site of action compared 
with using plasma concentration to drive the PD model. This 
study also adds to existing knowledge by providing a spe- 
cific application example validated against clinical data that 
confirms predictions of the discordant effect of transporter 
activity on the concentration of the statins in plasma and liver. 
It is anticipated that the approach used is applicable to other 
drugs with intracellular sites of action that rely on transporter- 
mediated processes for distribution to the site of action. It 
would be useful in the prediction of the impact of transporter- 
mediated drug-drug interactions on PD response in addition 
to the effect of genetic variations in transporter activity. 

METHODS 

Development of a PBPK/PD model describing the effect 
of rosuvastatin on cholesterol synthesis rate 

A PBPK/PD model of rosuvastatin in the north European 
Caucasian healthy volunteer population was constructed 
in the Simcyp Simulator (version 12 Release 2) as outlined 
in Figure 1. Detail of the PBPK model inputs for rosuvas- 
tatin has been described previously,^^ and further details of 
the calculation of the unbound concentration of a monopro- 
tic acid in the liver extracellular and intracellular water are 
given in the Supplementary Methods. Briefly, the dispo- 
sition of rosuvastatin was described using a whole-body 
PBPK model with tissue partition coefficients predicted by 
the method of Rodgers and Rowland,^** assuming perfusion- 
limited distribution for tissues other than the liver and gut. 
Distribution of rosuvastatin to the liver was described by 
a permeability-limited liver model that included the trans- 
porter-mediated intrinsic clearance (CL^,.^) for the sinusoidal 
uptake transporters 0ATP1B1, 0ATP1B3, 0ATP2B1, and 
sodium-dependent taurocholate cotransporting polypeptide 
and the canalicular efflux transporter breast cancer resis- 
tance protein. The model assumes that there is no trans- 
porter-mediated basolateral efflux from the liver, although 



recent work has indicated a role of basolateral efflux trans- 
porters hepatic efflux. Absorption of oral rosuvastatin was 
modeled using the Advanced Dissolution, Absorption and 
Metabolism (ADAM) modeP" and included active efflux by 
breast cancer resistance protein. 

Parameters describing rosuvastatin absorption, distribu- 
tion, metabolism, and elimination were not changed from 
the previously published values, with the exception of 
the 0ATP1B1 hepatic sinusoidal uptake transporter intrin- 
sic clearance for the three 0APT1B1 c.521T>C genotypes. 
Using published clinically observed concentration-time 
data,^^ the Simcyp Parameter Estimation module was used 
with the Nelder Mead optimization algorithm to obtain the 
uptake clearance of rosuvastatin into the liver for 0ATP1B1 
genotypes with C.521TT, TC, and CC sequence variations. 
To reduce the likelihood of overestimating interindividual 
variability in liver disposition within the different 0ATP1B1 
genotype groups, the variability in 0ATP1 B1 relative activity 
in the liver was adjusted to maintain the same overall coef- 
ficient of variation (CV) for the 0ATP1 B1 uptake CL.^,.^ in the 
northern European Caucasian population as when genotype 
was not considered (CV 69%). Equal variability in 0ATP1 B1 
activity was assumed for each genotype, and an adjusted 
CV of 44% was calculated using Eqs. 1 and 2. 



-Nxl 



N 



CV(%)=100x- 



(1) 



(2) 



where is the overall weighted standard deviation for all 
groups, n. is the fractional frequency of genotype / in the 
population, a. is the standard deviation of genotype group /, 
x, is the mean value of genotype group /, N is the total fre- 
quency of all genotypes in a population (N = 1), and \s 
the weighted mean for the population. The frequency of each 
genotype was calculated based on a weighted mean for the 
frequency of the c.521T>C SNP of 17% in the northern Euro- 
pean Caucasian population, and assuming Hardy- 
Weinberg equilibrium. 

The structural model for the effect of rosuvastatin on cho- 
lesterol synthesis was coded using the custom PD script- 
ing facility within the simulator (see Supplementary Data) 
and was based on the report by Aoyama et al.^^ Equations 
describing the MVA concentration in plasma were used as 
reported in this publication, with the exception of the con- 
centration input to the PD model (see Supplementary 
Methods). The parameters for baseline MVA concentration 
in plasma were kept as in the original publication. However, 
the drug effect (inhibition of MVA synthesis) in our model 
was driven by the unbound intracellular water concentration 
(CU|^,,) in liver, as opposed to plasma in the original model. 
Therefore, associated values (IC^^ and the Hill coefficient for 
the inhibitory sigmoid E^^^ function) were obtained by refitting 
the data using the Simcyp Parameter Estimation module with 
the Nelder-Mead optimization algorithm and clinical data for 
the change in MVA concentration for morning and evening 
dosing of rosuvastatin in dominantly wild-type 0ATP1B1 
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genotypes (using the compound file default 0ATP1B1 CL^,.^ 
109 ^il/min/10'* hepatocytes)^^ 

Although parameter variability to account for interindividual 
variability in plasma and tissue concentration profiles was 
incorporated in the PBPK model, this was not possible for the 
PD model because the published PD model adapted for this 
work was generated from fitting of average profiles. 

Simulations 

All simulations used the Simcyp north European Cauca- 
sian healthy volunteer population, and 10 trials were simu- 
lated. Where simulations aimed to replicate clinical studies, 
the trial design was matched as closely as possible to the 
study population in terms of dosing regimen, trial size, subject 
age range, and the proportion of females, as summarized in 
Supplementary Tabie S2. For predictive simulations, all sim- 
ulations used a lOmg oral dose of rosuvastatin dosed daily 
for 5 days and a trial size of 10 subjects with an age range of 
20-50 years and the proportion of females as 50%. A lOmg 
dose was selected because this was used by both the PK 
and PD study used in model development. ^'^^ The simulated 
plasma concentration and PD response reached steady state 
prior to the fifth dose of rosuvastatin (data not shown), thus 
the results for the final dose reflect predictions at steady state. 

The overall PD effect was summarized by calculating the 
% reduction in 24-h IVIVA AUC from baseline at steady state 
using the trapezium rule and as described previously.^^ 

Sensitivity analysis for the impact of total uptake 
transporter CL.^,^ on rosuvastatin disposition and PD 
response 

The automated sensitivity analysis option within the Simcyp 
Simulator was used to investigate the effect of total uptake 
transporter CL^^^ on the plasma, liver, and muscle exposure 
to rosuvastatin and the change in PD response. The specific 
output variables investigated were plasma, liver CU|„, and 
muscle AUC(,_24f, and % reduction in MVA AUC from baseline. 
Sensitivity analysis was selected to investigate overall CL^,.^ 
over the range of normal transporter activity for rosuvastatin 
uptake (250 pl/min/1 0" cells; approximately the sum of uptake 
transporter activity as reported in Jamei et a/.") to a com- 
plete loss of uptake transporter activity (0 |jl/min/10'' cells). 

The El is a dimensionless expression of sensitivity that 
measures the relative change in an output variable Q (e.g., 
AUC) for a relative change in the input parameter P (in this 
case, CL^^^)."^ El was calculated as follows for each input 
parameter value n: 



EL = 



OiPn] 



SI 



(3) 



where Q(P„) is the value of Q when P= P„ and SI is a mea- 
sure of the change in the output variable Q per unit change 
in the input parameter value P„ ^ , from its initial value P^, 
calculated using Eq. 4. 



Sl„ = 



Q(P„,0-Q(P.) 



(4) 



Pn + -\ Pn 

For the highest input parameter value, N, SI is approximated 
as follows: 



Sl/v - 2Slw-i -Sl/v_2 



(5) 
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Study Highlights 

WHAT IS THE CURRENT KNOWLEDGE ON THE 
TOPIC? 

^ Disposition of rosuvastatin to its target organ, the 
liver, involves the action of several uptake trans- 
porters, including the polymorphic transporter 
0ATP1B1. Reduced hepatic uptake of statins 
has been shown to increase plasma concentra- 
tion but has little effect on liver concentration. 

WHAT QUESTION DID THIS STUDY ADDRESS? 

^ The question addressed was whether using a 
PBPK model predicted unbound liver concen- 
tration as opposed to the plasma concentration 
to drive the PD model would improve predictions 
of the effect of genetic variation of 0ATP1B1 
activity on the therapeutic response. 

WHAT THIS STUDY ADDS TO OUR KNOWLEDGE 

This study demonstrates a method that allowed 
improved ability to predict the impact of the 
activity of transporters involved in uptake to the 
site of action on the PD effect. 

HOW THIS MIGHT CHANGE CLINICAL 
PHARMACOLOGY AND THERAPEUTICS 

i/ The approach used is applicable to other drugs 
that rely on transporter-mediated processes for 
distribution to an intracellular site of action and 
offers a means to better differentiate variability 
resulting from drug disposition from true PD 
variability. 
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